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We study damage spreading among the opinions of a system of agents, subjected to the dynamics 
of the Krause-Hegselmann consensus model. The damage consists in a sharp change of the opinion 
of one or more agents in the initial random opinion configuration, supposedly due to some external 
factors and/or events. This may help to understand for instance under which conditions special 
shocking events or targeted propaganda are able to influence the results of elections. For agents 
lying on the nodes of a Barabasi- Albert network, there is a damage spreading transition at a low 
value €d of the confidence bound parameter. Interestingly, we find as well that there is some critical 
value e s above which the initial perturbation manages to propagate to all other agents. 
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Most natural, social and technological systems are con- 
tinuously subjected to external stimulations of all kinds. 
Examples are infections, hacker or terrorism attacks, 
noise, errors, etc. It may happen that, due to these 
shocks, a sudden change occurs in some feature of a lim- 
ited number of subjects (thinking about human beings 
the feature could be health, religion or, like in this paper, 
opinion) , and that successively, through interactions with 
other members of the system, this perturbation spreads 
until it eventually affects a big portion of the system. The 
study of such processes is of great importance, in order 
to take the effects of eventual future local failures under 
control, but it is also an important method to investigate 
the dynamics of a system. 

Damage spreading (DS) was originally used by Kauff- 
man as a tool for studying biologically motivated dy- 
namical systems. In physics, the first investigations fo- 
cused on the Ising model 2] . Here one starts from some 
arbitrary configuration of spins and creates a replica by 
flipping one or more spins; after that one lets both con- 
figurations evolve towards equilibrium according to the 
chosen dynamics under the same thermal noise (i.e. iden- 
tical sequences of random numbers). It turns out that 
there is a temperature T^, near the Curie point, which 
separates a phase where the damage heals from a phase 
in which the perturbation extends to a finite fraction of 
the spins of the system. 

Meanwhile there is a sizeable literature on this prob- 
lem, which finds applications in many fields of modern 
science. Epidemiology, for example, is by definition the 
study of a particular DS problem [j]. After the recent 
discovery that many systems in nature and society can 
be described as complex networks, scale free graphs have 
been intensively investigated and many classical prob- 
lems have been reformulated on such special topologies 
0. In particular it is very interesting to understand the 
mechanisms by which diseases, information, computer 
viruses, etc. spread over networks. 

In this paper the network represents the system of ac- 
quaintances between people and we study the following 
problem: suppose we have a community of voters at the 
beginning of an electoral campaign, during which the vot- 



ers shape their own opinions through relationships with 
their friends. If a small set of voters for any reason sud- 
denly change their mind at this initial stage, would the 
final outcome of the election be influenced and, if yes, 
to which extent? Very recent history delivers a dramatic 
example: the shock caused to the Spanish people by the 
bombs in Madrid on March the 11th 2004 turned over the 
outcome of the national elections, which seemed already 
decided till that day. 

What we need is a model that describes how people 
convince each other. In the last years quite a few models 
of opinion dynamics have been proposed, like those of 
Galam Deffuant et al. Q, Krause-Hegselmann (KH) 
HJ], or Sznajd 0, and sociophysics simulations have be- 
come a fruitful field of research 0, [nj . As far as our 
DS problem is concerned, some results on the Sznajd 
model recently appeared [12], but for an improbable so- 
ciety where agents sit on the sites of a square lattice. 

Here we present the first systematic study on the sub- 
ject. Our agents are on the nodes of a Barabasi- Albert 
(BA) network [l^, which represents a more realistic 
model for the structure of social relationships. The net- 
work is constructed by means of a growth process start- 
ing from m nodes which are all connected to each other. 
Nodes are then added one by one and each of them forms 
m edges with the existing vertices, such that the proba- 
bility to get linked to a node is proportional to the num- 
ber of its neighbours. At the end, the number of agents 
with degree k, i.e. having k neighbours, is proportional 
to 1/k 3 for k large, independently of m. Our simulations 
show that the main results are only weakly dependent on 
m, so we focused on the case m = 3. 

We adopted the opinion dynamics of the KH consensus 
model. In this model, the opinions are real numbers be- 
tween and 1 and a confidence bound parameter e, also 
real in [0:1], is introduced. One starts from a random 
distribution of opinions. If we want to update the opin- 
ion Si of agent i we have to select among all neighbours of 
i only those agents whose opinions are compatible with 
Si, i.e. those agents j such that \si — sj\ < e; next, i 
takes the average opinion of the compatible agents. We 
chose to update sequentially the status of the agents, in 
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an ordered sweep over the whole population; in this way, 
the dynamics does not require random numbers and is 
therefore truly deterministic, at variance with other con- 
sensus models. At some stage, the system will converge 
to a configuration that the dynamics is unable to modify. 
This configuration represents the equilibrium state and 
the final opinion distribution, which is given by a set of 
(5-functions, depends on e. On BA networks, we found 
that the threshold for complete consensus is e c = 1/2, 
independently of m. 

For the DS analysis we followed the procedure that 
we exposed above for the Ising model. After creating a 
random opinion configuration, we produced a replica of 
it in that we changed the opinion of a single agent i. The 
results do not depend on the exact number of perturbed 
sites, as long as they are a fraction of the population N 
that vanishes in the limit when N goes to infinity |14| . 
We perturbed the opinion as follows: if Si > 1/2 the new 
opinion becomes Sj — 1/2, otherwise Sj + 1/2. 



(the time unit is one sweep over the network). There is 
a characteristic pattern with an initial phase in which 'F 
grows rapidly, followed by a very slow relaxation to the 
final value. 
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FIG. 1: Time evolution of the fraction of perturbed opinions 
for several values of e. Here the population is N = 10000. 

After each sweep over the network we calculated the 
fraction 'F of different opinions in the two configurations. 
The simulation stops when both systems reach their final 
states, which happens when no agent changed opinion 
during an iteration. In all our calculations we used double 
precision real numbers and we decided that two opinions 
are the same if they differ from each other by less than 
10~ 9 , otherwise they are different. For the results to 
have statistical significance we took averages over 1000 
samples for all values of e and of the population N. 

We stress that on our network, like in a real society, 
some nodes are more important than others. If we aim 
at spreading an opinion in the society, we should better 
try to convince people with many friends than persons 
with few social contacts. We then expect that the per- 
turbation on the whole system will be more relevant if 
we initially damage a hub than a loosely connected node 
and we investigated both situations. 

First, we analyzed the case in which the perturbed 
node is a hub. Fig. ^ shows the variation with time of *F 



FIG. 2: Variation of \& with the confidence bound parameter 
e. 

As we can see from the figure, the relaxation time in- 
creases up to e ~ 0.1, then it decreases. At the end of 
relaxation the number of damaged sites stays constant 
although the system keeps evolving. This is due to the 
fact that, at advanced stages of the evolution, the opin- 
ions are grouped in clusters. Agents of different clusters 
cannot interact with each other, as their opinions differ 
by more than e. Because of that, if a cluster contains 
perturbed opinions, most of its agents will be sooner or 
later affected [lfij . otherwise it can never be reached by 
the perturbation. 

In Fig. |2 we plot f as a function of e, for different 
network sizes. We see that the damage rises fast with e 
and that, for e larger than about 0.05, more than half 
of the agents have been affected. The inflection of the 
curves at low values of the confidence bound hints to the 
existence of a DS transition, like in the Ising model. If 
there is indeed a phase in which the damage affects only 
a vanishing fraction of the system, we should find that 
there 'F goes to zero when TV goes to infinity. We then 
looked for scaling behaviour of 'F with N. In Fig. [31 
we plot as a function of N for several values of the 
confidence bound. We can see that the points can be 
quite well fitted by a simple power law up to e ~ 0.015. 
For higher values, a saturation to a non-zero \& takes 
place. As estimate of the DS threshold e<j, we took the 
value which gave the smallest \ 2 f° r the fit with the power 
law: we found e e i — 0.013(3). The error marks the range 
where the above-mentioned x 2 (P er degree of freedom) is 
below 1. 

Coming back to Fig. |3 we notice that in contrast to 
the Ising or Kauffman models, we have a second thresh- 
old e s above which all agents are affected by the initial 
local perturbation. To our knowledge, this is a truly 
novel feature for a DS process, and we expect it to hold 
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FIG. 3: Dependence of ^ on the population N for low values 
of e. 



as well for the opinion dynamics of Deffuant et al. By 
studying the dependence of this threshold on N, we ex- 
trapolated its infinite- AT limit e s = 0.500(1). We remark 
that this is just the threshold for complete consensus of 
our system. So, if we are in the consensus regime and 
a single agent suddenly changes its mind, this suffices to 
(slightly) modify the final dominant opinion of the total 
population. 

We tried to check whether this unexpected feature is 
specific of the particular social topology we have chosen, 
or whether it is exclusively due to the dynamics. Simula- 
tions of a society where agents sit on the sites of a square 
lattice, with periodic boundary conditions, confirm the 
result and e s is again 1/2. We have examined as well a 
community where each agent has relationships with all 
others; here we have that e s is about 0.07, much lower 
than the consensus threshold 0.21. For this special so- 
ciety there is no DS transition, as each agent effectively 
interacts with a finite fraction of the system, and it can 
be proved that £d(N) vanishes when N — > oo. 

We remind that we have introduced the damage right 
at the beginning of the evolution. We have also per- 
formed some tests to check what happens if we instead 
perturb the system after some evolution steps. Now 
the perturbed agent stays in a society where people are 
mostly divided in groups of close opinions, and such com- 
munities will evolve separately from each other. In this 
way the shocked agent can interact only with a smaller 
portion of the system, i.e. with a few clusters of agents, 
and the amount of damage will drop. On the other hand, 
we find that both and e s remain the same. 

Furthermore we have studied the effective opinion vari- 
ations of the agents if damage is introduced. For this 
purpose we divided the opinion of each agent in the per- 
turbed configuration by the corresponding value in the 
unperturbed configuration. In Fig. 0]we show the distri- 
bution of the opinion ratios at the end of the time evolu- 
tion, for two values of the confidence bound. As we can 
see, when e is small, the distribution is strongly peaked 
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FIG. 4: Histogram of the ratios of the agents' opinions with 
and without damage. The confidence bound is e = 0.05 (top) 
and 0.45 (bottom); the population is N = 1000. 



at the value 1, and the opinions vary continuously though 
very little in most cases. As we approach the consensus 
threshold, instead, we notice that the opinion variations 
are no longer continuously distributed, and other narrow 
peaks appear, which shows that most values of the ratios 
are suppressed and discontinuous jumps, corresponding 
to drastic opinion changes, are allowed. 

Let us now check what happens when the initially dam- 
aged agent sits on a node with low degree. In Fig. 
we compare the DS curve obtained in this case with the 
curve for a perturbed hub. We see that, for any e, $ is 
larger when we damage the hub, as expected. On the 
other hand we find that both and e s are the same 
as before. That relies on the small world effect on 
scale free networks like ours. In fact, each node can 
reach any other through a small number of intermedi- 
aries. In this way, even if we perturb a loosely connected 
node, within the first evolution steps the perturbation 
will have reached quite a few nodes with much higher 
degrees, which brings us back to the previous case. The 
damage is larger if we perturb a hub because more agents 
can be reached at the beginning of the evolution; soon af- 
ter that, as we said above, clusters of opinions are formed 
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which do not interact with each other and the perturba- 
tion can exclusively spread within the affected clusters. 
In conclusion, due to the small world effect, \& and the 
fraction $ = I — 'F of unperturbed agents are of the same 
order of magnitude in both situations: if 'F ($) is zero in 
one case, it will be zero in the other case too. 
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FIG. 5: Comparison of the DS curves corresponding to the 
initial perturbation of a node with low degree (squares) or a 
hub (circles). The number of agents is N = 1000. 

We have studied damage spreading for the Krause- 
Hegselmann opinion dynamics on Barabasi- Albert net- 
works. We distinguish three phases in the confidence 



bound space, corresponding to zero, partial and total 
damage, respectively. The existence of a phase where the 
perturbation affects all agents is new for damage spread- 
ing processes and is independent of the social topology, 
so it only relies on the dynamics. This feature seems 
unrealistic, but we should consider that our dynamics is 
very simple [l7| . Moreover, we let our system reach in 
any case its final state, but the evolution time grows with 
the size of the community and would be very long for a 
realistic population of voters; normal electoral campaigns 
last a few months, so we should interrupt the evolution 
process at an earlier stage and the damage could then be 
limited. The probability to have large variations of the fi- 
nal agents' opinions induced by the initial perturbation is 
rather small, but it increases with e and, by approaching 
the consensus threshold, forbidden bands appear. The 
amount of the damage depends on the degree of the dam- 
aged node but the thresholds for damage spreading and 
saturation do not, because of the small world effect. 

In the future it would be helpful to use other opin- 
ion dynamics to check for the consistency of the results. 
For more realistic analyses of the problem, together with 
eventual refinements of the existing consensus models 
and of the social topology, it is important to include as 
well other factors in the dynamics, like advertising and 
noise. 
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